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We compute the ground-state properties of finite systems of neutrons in an external harmonic 
trap, interacting via the Minnesota potential, using the "exact-exchange" form of orbital-dependent 
density functional theory. We compare our results with Hartree-Fock calculations and find very 
close agreement. Within the context of the interaction studied, we conclude that this simple orbital- 
dependent functional brings conventional nuclear density functional theory to the level of Hartree- 
Fock in an ab initio fashion. Our work is a first step towards higher-order ab initio nuclear functionals 
based on realistic nucleon-nucleon interactions. 



I. INTRODUCTION 

Density functional theory (DFT) is a general theory of 
quantum many-body systems with a long history. In its 
modern version, that began with the work of Hohenberg, 
Kohn and Sham [1, 2], DFT has become an essential tool 
in quantum chemistry as well as in materials science and 
condensed matter physics (sec e.g. [3-5]). This is largely 
due to the fact that large systems, intractable by ab initio 
methods like Coupled Cluster (CC) or Quantum Monte 
Carlo (QMC) even with modern computational power, 
are typically within the reach of Kohn-Sham (KS) DFT 
with relatively modest computational resources [6]. In 
spite of significant progress and multiple successes, the 
fundamental challenge remains the same in all applica- 
tions: the central object of the theory, namely the energy 
density functional (EDF) is a priori unknown and must 
be built either phenomenologically or from first prin- 
ciples (or some combination thereof) by implementing 
some kind of approximation scheme. Conventional DFT 
has slowly evolved from the first simple and often uncon- 
trolled approximation strategies of the early days, such 
as Thomas-Fermi theory [7] and the Local Density Ap- 
proximation (LDA) [2], to more sophisticated semi-local 
approaches including the Generalized Gradient Approxi- 
mation (GGA), meta-GGAs and hyper-GGAs [8]. 

Progress in the field of quantum chemistry [9] has 
pushed the boundaries of such conventional approaches 
to DFT by allowing electronic EDFs to depend explicitly 
on the single-particle orbitals of KS DFT. While the set 
of possible functionals is thus enlarged, this additional 
freedom comes at the price of increased formal complex- 
ity and computational demand. Indeed, the determina- 
tion of the KS auxiliary one-body potential v KS in the 
case of orbital-dependent functionals necessitates the so- 
lution of the Optimized Effective Potential (OEP) inte- 
gral equation [10, 11], which prevented widespread use 
of these functionals for a long time. Significant advances 
in the last decade, however, have enabled the numeri- 
cal solution of the OEP equation in a more straightfor- 
ward and systematic fashion, thereby allowing orbital- 



dependent functionals to enter the mainstream of elec- 
tronic DFT [12, 13]. 

These developments were often motivated by the in- 
adequacies of conventional electronic EDFs concerning 
some practical issues (such as their inability to deal 
with Van der Waals forces, to predict the existence of 
negative ions and to properly account for electron self- 
interaction [9]) as well as formal problems (such as the 
lack of particle-number derivative discontinuities [14], in- 
ability to reproduce the correct long-range tail of the one- 
body potentials, etc. [9]). Orbital-dependent functionals 
have not only mitigated these issues, all of which stem 
from a poor description of particle exchange in GGA- 
based EDFs, but they have also cleared a path towards ab 
initio DFT. Indeed, many-body perturbation theory ex- 
pressions for the total energy (at first order in its simplest 
form, i.e. Hartree-Fock, but also at second and higher 
orders, and in resummed forms) are generally orbital- 
dependent functionals, which have a clear and direct con- 
nection to the microscopic Hamiltonian [15]. 

Many-body systems of electrons in external fields are 
the standard arena for the application of DFT. However, 
nuclear physics shares the same interests in many-body 
techniques and faces very similar challenges. DFT has 
therefore become a standard tool for the computation 
of the properties of heavy nuclei. In nuclear DFT, the 
electrons are replaced with nucleons, and the nucleon- 
nucleon interaction plays the role of the Coulomb inter- 
action. The vast majority of nuclear functionals, whether 
phenomenologically motivated or derived from first prin- 
ciples, are largely in the category of GGAs (broadly de- 
fined). The parameters in these functionals are usually fit 
to a subset of stable nuclei (with the exception of some 
of the latest functionals, such as SLy4 or SkO, which 
have been fit to some experimentally accessible unsta- 
ble nuclei) and are therefore not directly related to our 
understanding of the nucleon-nucleon interaction. 

The main objective of this work is to explore the 
prospects of ab initio DFT for the nuclear case. A more 
complete survey of the current status of the field can 
be found in Ref. [16]. In this work, we implement one 
of the simplest possible orbital-dependent DFTs derived 
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from an underlying hamiltonian, the so-called "exact- 
exchange" (EXX) form. We apply this to the case of neu- 
tron drops interacting via the Minnesota potential [17]. 
We then compare the results of this DFT with exact 
Hartree-Fock (HF) for various numbers of neutrons in 
a harmonic external potential. 

The EXX functional is easily defined: it consists of the 
HF energy, in which the single-particle HF orbitals are 
replaced with the KS orbitals. Thus, the fundamental 
difference between EXX DFT and exact HF is that at 
each step in the iterative optimization procedure the HF 
approach involves a non-local auxiliary potential, whereas 
KS DFT (regardless of the form of the functional) uti- 
lizes a local auxiliary potential v KS . (It remains an open 
question whether such local auxiliary field will be too 
constraining for nuclear physics, see e.g. Ref. [18].) In 
this sense, EXX DFT can be regarded as a constrained 
optimization of the HF energy, where the constraint con- 
sists in demanding the locality of v KS . As a consequence, 
the ground-state energies obtained via EXX DFT should 
be expected to be higher than those of HF. As we shall 
see, the differences are very small in the case we study, 
but this will in general depend on the form of the inter- 
action. 

While EXX DFT is formally simple, we do not claim 
that it is accurate in an absolute sense, but rather that 
it represents an extremely accurate approximation to 
exact HF. As mentioned above, this property will in 
general depend on the interaction. However, we wish 
to stress this point as a promising feature of orbital- 
dependent DFT, even though HF calculations are well 
known to be a poor description of nuclei due to the non- 
perturbative nature of the nuclear interaction at short 
distances (c.f. high momenta). A more realistic descrip- 
tion must both transcend perturbation theory and in- 
clude pairing correlations at least at the mean-field level, 
i.e. a la Hartree-Fock-Bogoliubov (HFB). Work in this 
direction has recently involved efforts towards taming the 
nuclear interaction at high momenta using Renormaliza- 
tion Group (RG) transformations [19] [24]. These are 
transformations that leave observables unchanged but 
reduce the strength of the potential at large momenta. 
Such transformations render the problem more pcrturba- 
tive while maintaining the hierarchy of many-body forces. 
In this context, our work may be regarded as a neces- 
sary first step in an ab initio DFT program that connects 
microscopic RG-transformed Hamiltonians with orbital- 
dependent DFT based on HFB plus perturbation theory. 



II. KOHN-SHAM DFT AND THE OPTIMIZED 
EFFECTIVE POTENTIAL 

In quantum chemistry, the "optimized potential 
method", or simply the "optimized effective potential", 
refers collectively to the use of orbital-dependent EDFs 
and to the determination of the KS auxiliary potential 
^ks by solving the OEP integral equation. For complete- 



ness, and in order to set our notation, we present here 
a short derivation of this equation, along with a brief 
review of KS DFT. For simplicity, we shall restrict our- 
selves to functionals that do not depend on the KS eigen- 
values; the corresponding generalization is easy to carry 
out. 



A. Derivation of the OEP equation 

The central tenet of DFT is the Hohenberg-Kohn (HK) 
theorem, whereby the existence of an energy density func- 
tional is asserted, of the form 

E[p]=F[p}+E eKt [p}, (1) 

where 

E cxt [p] = J dx u ext (x)p(x), (2) 

such that F[p] depends only on the one-body density 
p(x). One may think of the latter as the total density, but 
in general it may denote spin, isospin, kinetic or anoma- 
lous densities, in which case x represents a collective in- 
dex for the coordinate and every other degree of freedom. 

The external potential w cxt represents the electric field 
of the ions, and it confines the system to a particular 
spatial region. This is a fundamental difference between 
the electronic and nuclear cases, since in the latter the 
system is self-bound, i.e. there is no v cxt [16]. 

According to the HK theorem, the functional form of 
F is determined solely by the interactions and not by the 
external potential w cxt ; the functional F is therefore said 
to be universal. The HK theorem is an existence theorem 
and gives therefore no instructions as how to build or find 
this functional. Kohn-Sham DFT takes a first formal step 
towards the explicit construction of F\p] by separating it 
into a non-interacting piece, i.e. the kinetic energy of the 
free system, and everything else: 

F[p]=T s + E in M (3) 

where 

^=EE/ dx *ux) (-^) *u*). (4) 

It should be stressed that T s represents the kinetic energy 
of the auxiliary KS system, which in general is different 
from that of the many-body system. Their difference is 
assumed to be accounted for in E int [p]. 

The ip kt7 (x) are a set of auxiliary single-particle orbitals 
(the KS orbitals), such that p = p^ + p^, where 

p CT (x) = f> feCT (x)| 2 , (5) 
fe=i 
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and N a is the particle number for spin a. Kohn-Sham 
DFT then proceeds to optimize the energy functional by 
solving a Schrodinger-like equation: 



fi 2 V 2 
2m 



+ «KS^( X ) 



VfcrW = CfcVfcaW. (6) 



for all k. As mentioned above, this should be contrasted 
with the HF approximation, in which the correspond- 
ing Schrodinger equation involves a non-local potential. 
In this sense, EXX DFT can be regarded as resulting 
from an HF minimization procedure, with the added con- 
straint that the auxiliary potential be local. In spite of 
this constraint, orbital-dependent DFT beyond EXX has 
the potential to surpass HF, and in fact does so in prac- 
tice in the electronic case (see e.g. Ref. [9]). 

By definition, the KS potential in Eq. (6) is given by 



where 



SV 



V = E n 



v ext ( x ) + c ( \ ' 



(7) 



(8) 



The second term in the r.h.s. of Eq. (7) will in general 
depend on the KS orbitals, such that Eq. (6) is to be 
solved self-consistently by starting with a guess for the 
orbitals or for v KS . 

It is at this point that orbital-dependent DFT departs 
from GGA DFT, in which functionals depend explicitly 
on the density and its gradients. Indeed, once we allow 
E int to depend explicitly on the KS orbitals, it becomes 
unclear how to determine the KS potential using its defi- 
nition Eq. (7). One of the simplest ways to proceed is to 
use the chain rule of functional differentiation and con- 
sider the following identity 



SV 



5V 5<p ktr ,(y) 



EE* 



+ c.c. 



(9) 



In order to proceed we vary both sides of the eigenvalue 
equation Eq. (6) (and its complex conjugate): 



Gr(y,x)^ CT (x), 



where 



with 



Gr(y,x) = ^G b (y,x), 



(10) 



(11) 



q^k ek e i 



is the Green's function. Equation (9) then becomes 



7Z\ = 1^ / dy Kn ^ G fc^(y> x )^fca( x ) +c-c. 



Sv 



(13) 



On the other hand, the HK theorem allows us to take 
implicit derivatives with respect to the density p CT (x): 



6V 



5v ksA x ) 



= / dx! 



fy<r( x l) ^KS,rr( X ) 



(14) 



( X l) 



5/0,7 (Xl) 



E / rfx l W KS,a( x l) ( PL( x l)G fcCT (xi,x)( / 9 fea (x) 



fc=l 



+ C.C, 



where we have used Eqs. (5) and (7). Combining 
Eqs. (13) and (14) we arrive at 



E (CWfeW + c - c -) = °. 



(15) 



*;=i 



where we have defined the "orbital shift" ip krT by 

SV 



V>L( X ) = / dxi 



/ 



x G feCT (xi,x). 



^KS,a( X l)^L( x l) 



(16) 



Equation (15) is the OEP integral equation that defines 
the KS potential of orbital-based DFT. 



B. Solving the OEP equation 

While there are in principle multiple ways to solve the 
OEP equation, we have found the one originally due to 
Kiimmel and Perdew [13] to be particularly useful. We 
shall refer to this method as the Kiimmel-Perdew (KP) 
algorithm. This algorithm is iterative and starts by not- 
ing that the orbital shifts can be found by solving a dif- 
ferential equation, namely 



fi 2 V 2 
2m 



+ «KS, ff ( x ) - £ fc 



and 



CW = A fc CT ( x VL( x ) 

(17) 



-J dx l(«KS^( X l)- U fea( x l)) ( f^L( x l)r'fc CT ( x l)(18) 

with 

1 SV 



</>L( x ) ^ fcCT (x) : 



(19) 



where we assume some form for the initial v KS a (x) . Once 
the orbital shifts are determined, v KS a is updated ac- 
cording to 



^ e s w CT (x)=< 1 s d CT (x) + c5 (7 (x), 



(20) 
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where c is a positive constant and the function 



and 



S ff (x)E^(C(xK„(x)+c.c.) (21) 



fc=i 



(c.f. Eq.(15)) is used as a local measure of the deviation 
of the current v KS with respect to the true solution of the 
OEP equation. Indeed, wherever t> KS(T (x) differs from 
the true OEP, S a (x) will be non-zero. The updated v KS 
can be used to re-compute the orbital shifts, keeping the 
KS orbitals and eigenvalues constant, until the desired 
convergence criterion is satisfied. The final i> KS is then 
re-inserted in the KS equation (6) to compute a new set 
of orbitals, closing the self-consistency loop. 

It is easy to see that Eq. (17) is singular. Adding 
any multiple of ip ka (x) to ip ka {x) gives a new non-trivial 
solution. As first noted in Ref. [13], however, Eq. (17) 
can still be solved in the subspace of interest by using 
the method of Conjugate Gradients [21]. Further details 
of this approach can be found in Ref. [13]. 

The determination of the orbital shifts ^* ka via Eq. (17) 
represents an additional computational cost that is not 
present in LDA- or GGA-based DFT. Indeed, in the lat- 
ter v KS rT is determined simply by inserting the density 
and its gradients in a pre-calculated analytic expression 
derived from the EDF via Eq. (7). In contrast, in the 
OEP method (within the KP approach) the calculation 
of all the V'fecr's necessitates the solution of N differential 
equations, where N is the number of particles. While 
this is a considerable increase in computational demand, 
it should be kept in mind that the ip^'s are independent 
from each other, such that the N differential equations 
can be solved in a parallel fashion with perfect scaling (up 
to communication costs required to broadcast % s (x) at 
the beginning and to compute 5 CT (x) at the end). 



III. EXX ENERGY DENSITY FUNCTIONAL 
AND THE MINNESOTA POTENTIAL 

We have used in this work the Minnesota nucleon- 
nucleon interaction of Ref. [17], which is given for pure 
neutron systems by 



V(r) = (V R (r) + V s (r)V s + V t (r)V t ) 



(22) 



where V s (Vt) is the operator that projects on the spin 
singlet (triplet) state and P r is the coordinate exchange 
operator. The potentials in the various channels have 
Gaussian forms given by 



V R 

v s 



V r oR,exp(-K R r 2 ) 
-Vbt exp(-K t r 2 ) 
-V 0s exp(-K s r 2 ), 



V 0R = 200.0 MeV, 
Vbt = 178.0 MeV, 
V 0s = 91.85 MeV, 



k r = 1.487 fm" 
0.639 fm" 



: 0.465 fm" 



(24) 



Our main reason for using this potential is that it pro- 
vides a semi-realistic yet easy-to-implement interaction. 
In addition, this potential is moderately soft and as such 
it is comparable to those obtained using renormalization- 
group methods [19]. Other approaches to ah initio 
DFT are also being explored using this simple potential, 
and No-Core-Full-Configuration (NCFC) and Coupled- 
Cluster (CC) results have recently become available as 
well, which makes this a good test case for OEP meth- 
ods. 

The energy density functional is of the form given in 
Eq. (3) and the mass in the kinetic energy term (Eq. (3)) 
has been set to m n = 939 MeV. In practice, one combines 
the kinetic term with the external Harmonic Oscillator 
(HO) potential energy. This is convenient because the 
KS orbitals can be expanded using the HO basis {4>j}i 



3 = 1 



a jk ,a = J dx <^*( X Vfc CT ( X )> 



(25) 
(26) 



and in this basis the sum T s +E cxt takes a simple diagonal 
form. 

In the EXX case considered in this work the interaction 
enters through 



(27) 



ijkl 



where we have used collective indices i to denote the pair 
(i,(Ti) for the basis element and the spin. Here, 



Pij = Picrfjaj — 5<7 z <jj ^2 a ik,a z a 3 k,a z ( 2 ^) 



fe=l 



is the one-body density matrix, and 



v m = v m -v m 



(29) 



(23) 



are the anti-symmetrized matrix elements of the inter- 
action. As mentioned above, this form of F[p] is simply 
that of the HF energy, but with the KS orbitals replacing 
the HF single-particle wavefunctions. The interaction in 
Eq. (23) is local and we therefore obtain 

Vijki= J rf x i rfx 2^( x i)^( x 2) y ( r )^( x i)^i( x 2)' ( 30 ) 

where r = \x 1 — x 2 | and the total spin of the initial (ik) 
and final (jl) states is the same. 
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IV. RESULTS AND CONCLUSIONS 

We have considered systems of 8 and 20 neutrons for 
various trap frequencies and basis sizes, and computed 
the total energy 



^kin ~t~ ^cxt T J^jj ' -^Fj 



where E kin = T s , and 



f/ cxt = J dx u ext (x)p(x), 



where 



2„2 



«cxt( x ) = 2 m "^ X 



(31) 



(32) 



(33) 



is the HO external potential. 

The total energy can also be obtained from eigenvalue 
sum rules, which we have verified numerically, and are 
given in the HF case by 



E,. 



J^ffc - E H - E F , 



(34) 



where each term is evaluated using HF orbitals, and in 
the EXX-DFT case by 

E to t =^2ek+E H + E F -^2 f rfx %Siff (x)^(x), (35) 

fc,<T a J 

where the EXX-DFT orbitals should be used. Notice 
that in both sum rules the eigenvalue sums go over both 
spins. The Hartree and Fock energies are respectively 
given by 



^ H ~~ 9 51 VijklPkiPlj 



ijkl 



E F = 



lyv. 



2 / ijlkPkiPlj- 
ijkl 



(36) 



(37) 



We have also computed the r.m.s. radius -J (r 2 ), where 



(r 2 ) = jf[dx r 2 p(r) 
as well as the form factor, given by 
F(q) — 47T dr r 2 p(r)'- 



qr 



(38) 



(39) 



of convergence of HF and EXX DFT, with respect to the 
size of the basis N majL , by showing data for iV max = 64, 
125 and 216. Figures 2 and 3 show the density profiles 
and the corresponding form factors for a fixed external 
harmonic potential of frequency Ml = 10 MeV and basis 
size iV max = 27, 64, 125 and 216. Throughout this work 
we have chosen the value of Ml for the basis equal to 
that of the trapping potential. As seen in the plots, in 
all cases the convergence pattern of HF and EXX DFT 
as a function of AL,_„ is the same. At fixed iV m „, on 
the other hand, the convergence patterns as a function 
of the respective HF and KS iterations are significantly 
different from each other. We provide a partial summary 
of our results in Table I. 

We have implemented a parallel code (using OpenMP) 
for both the HF and OEP calculations. The latter use 
some of the HF routines and a set of extra ones to solve 
the OEP equations. Using 8 processors and some limited 
optimizations of the parameters (such as the constant c 
in the KP algorithm and the number of iterations in the 
KP loop) we have found that, at fixed accuracy for the 
convergence criterion, OEP calculations require about a 
factor of 2 — 3 more iterations than HF, and each OEP 
iteration takes about a factor of 2 — 3 more time than its 
HF counterpart, for the basis sizes and particle numbers 
studied here. Thus, overall the OEP method at the EXX- 
DFT level is a factor of 4 — 9 slower than HF. However, 
further optimizations are possible if for example MPI is 
used in the parallelization of the calculation of the orbital 
shifts. This would have significant impact in calculations 
for larger particle number and larger basis size, as well 
as for more sophisticated functionals. 
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where q = |q| and r = |x|. 

Our results for the internal energy E tot — {J ext per par- 
ticle, as a function of the r.m.s. radius, are shown in 
Fig. 1, for various frequencies of the trapping potential. 
Also shown in Fig. 1 are the results corresponding to 8 
and 20 neutrons, and in all cases we display the degree 



FIG. 1: (Color online) Internal energy as a function of the 
r.m.s. radius for 8 (lower curve) and 20 (upper curve) neu- 
trons, computed with HF and EXX DFT, for two different 
basis sizes, namely iV max =125 and 216. From left to right, 
the trapping potential corresponds to hQ = 20, 15, 10, 5 and 
3 MeV. 
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FIG. 2: (Color online) Density profiles of a system of 8 neu- 
trons in a HSl = 10 MeV trap as a function of the basis size 
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FIG. 3: (Color online) Form factor of a system of 8 neutrons 
in a HQ — 10 MeV trap as a function of the basis size -/V max , 
computed with HF and EXX DFT. 



As mentioned in the introduction, this work represents 
the first step in a program whose objective is to construct 
a nuclear energy density functional from first principles, 
employing orbital-based methods. Ideally, this functional 
would be able to at least qualitatively predict the prop- 
erties of heavy nuclei and would only be based on our mi- 
croscopic understanding of the internucleon interaction. 
We do not exclude the possibility, however, that an ac- 
curate quantitative description of energies and radii may 
require readjusting certain parameters in this functional. 
In this sense, our work is a step forward from LDA- and 
GGA-based DFT, which are typically phenomenological 
in nature in the nuclear case, and therefore limited in 
their ability to predict the properties of unknown sys- 
tems. Much work remains to be done, however, to bring 
this approach to a level of accuracy that is competitive 
with state-of-the-art nuclear EDFs. 

Orbital-based DFT is one possible road out among sev- 
eral currently pursued in the nuclear DFT community. 
An alternative way is to find a controlled scheme that 
approximates the non-local ab initio functional with a 
local one and thereby facilitates the application of stan- 
dard KS DFT (see e.g. [22]). 

The next step in this ab initio DFT program will be 
to treat nuclei by implementing realistic nuclear inter- 
actions such as chiral interactions softened with mod- 
ern renormalization group methods. This will require 
extending the current framework to include many-body 
forces and non-central as well as non-local terms in the 
interaction. This is currently underway; we have suc- 
cessfully performed proof-of-principle calculations with 
three-body forces, and shown that the formalism for 
many-body forces at the EXX level is a straightforward 
extension of our presentation in Sec. II [23]. Extensions of 
the OEP method currently in progress involve improved 
functionals, e.g. those generated by second-order per- 
turbation theory, as well as modifications to account for 



V. SUMMARY AND OUTLOOK 

In this work we have computed the energy, radius, den- 
sity profile and form factor of finite systems of neutrons 
in a harmonic trap. For this purpose we have used pure 
HF and orbital-based DFT in the EXX form, implement- 
ing the algorithm of Kiimmel and Perdew to solve the 
OEP equation. Our results for the OEP approach agree 
at the 0.1% level or better with those of pure HF for 
all the quantities we computed. We note in particular 
that such an agreement was reached for each basis size 
even though HF and OEP converge to their respective 
results in completely different ways. Although limited 
in scope, our work shows that it is possible to capture 
the exchange aspects of the nuclear interaction in an ex- 
tremely accurate fashion completely within the context 
of local KS-DFT. 



TABLE I: Summary of results for the energies in MeV for 
Hartree-Fock (HF) and Exact-Exchange DFT (EXX), broken 
up into total (E tot ), kinetic (E Mn ), Hartree {E H ), Fock (E F ), 
and internal (Z? t = E tot — £/ cxt ). Also shown is the r.m.s. 
radius \J (r 2 ). 



Method 


N 


hQ, 


Etot 




E H 


E F 


Ei 




HF 


8 


20 


296.36 


200.92 


-225.18 


159.35 


135.09 


2.0438 


EXX 


8 


20 


296.36 


200.99 


-225.24 


159.39 


135.14 


2.0435 


HF 


8 


10 


142.43 


111.22 


-125.38 


83.70 


69.53 


2.7484 


EXX 


8 


10 


142.44 


111.27 


-125.44 


83.74 


69.56 


2.7478 


HF 


8 


3 


44.51 


35.18 


-30.51 


19.05 


23.72 


4.8927 


EXX 


8 


3 


44.51 


35.19 


-30.52 


19.05 


23.72 


4.8921 


HF 


20 


20 


941.75 


707.69 


-920.03 


645.08 


432.74 


2.2965 


EXX 


20 


20 


941.77 


707.83 


-920.19 


645.20 


432.84 


2.2963 


HF 


20 


10 


456.25 


382.73 


-488.61 


326.55 


220.68 


3.1246 


EXX 


20 


10 


456.26 


382.88 


-488.80 


326.69 


220.77 


3.1240 


HF 


20 


3 


143.52 


119.55 


-118.19 


74.15 


75.52 


5.5958 


EXX 


20 


3 


143.52 


119.57 


-118.21 


74.17 


75.53 


5.5954 
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pairing correlations. 
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